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We show how spatiotemporal fluctuations can induce spontaneous symmetry breaking in systems 
which are perfectly symmetric in the absence of fluctuations. We illustrate this in the context 
of the autocatalytic production of chiral enantiomers from achiral reactants in reaction-diffusion 
systems. The mean field steady state is chiral symmetric; spatiotemporal fluctuations induce a 
novel (molecular) chiral ordering and sharp phase transitions including reentrance. We discuss its 
implications in the context of the emergence of molecular homochirality. 
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Fluctuations are known to change the nature of phase 
transitions or even destroy long-range order, and so it 
comes as a surprise that fluctuations may induce sym- 
metry breaking transitions in systems where such sym- 
metry breaking would not have occured in its absence 
Q . In this Letter we study the autocatalytic production 
of chiral enantiomers from achiral reactants in reaction- 
diffusion systems whose mean field steady state is chiral 
symmetric; spatiotemporal fluctuations however induce 
a novel chiral ordering and sharp phase transitions. 

While arriving at these results we comment on the con- 
ventional understanding of homochiral amplification in 
autocatalytic systems. F.C. Frank Q proposed a sim- 
ple chemical kinetics model for the amplified production 
of a chiral enantiomer from a racemic mixture seeded 
with an initial bias, based on autocatalysis of the favoured 
enantiomer and antagonism of the disfavoured. Follow- 
ing recent successful realisations of this scheme Q, it is 
widely believed that both an initial chiral seed and an- 
tagonism are necessary conditions for the amplification 
of homochirality in autocatalytic systems @ . Contrary 
to this, we show that complete chirality discrimination 
and exponential amplification can occur, (i) without an 
initial chiral seed and (ii) without specific antagonism, 
simply via stochasticity inherent in autocatalytic chemi- 
cal reactions. 

Consider a system of achiral molecules A and B which 
diffuse and react on collision to form chiral enantiomers 
R and R*, 



A + B ^ R 

fc 2 

A + B ^ R* 

*=2 



(1) 



The forward reactions are autocatalytic with rates k\ (R) 
and kx(R*) - existing enantiomers favour their own pro- 
duction - while the backward reactions with identical 
rates k 2 are non-autocatalytic. A simple linear model 
of autocatalysis is k\(R) = kiP(N R ), where P(N R ) is 
the probability of encountering R in the forward reac- 
tion. In our symmetric autocatalytic model, we take the 
constants k\ to be the same for the two reactions. The 



reverse reaction, on the other hand, has no reason to be 
autocatalytic; even so our qualitative conclusions remain 
the same were the reverse reactions autocatalytic too. 
Note that the molecules interact only through the auto- 
catalytic chemical reactions, which occur on collision. 

First define the mean-field model for this autocat- 
alytic reaction-diffusion system, valid when the number 
of molecules per unit volume is large so as to have a well- 
defined concentration field which varies smoothly and dc- 
terministically with space and time. The mean-field con- 
dition implies that local fluctuations in the concentration 
are much smaller than the mean, \/ n 2 — n 2 <C n. Under 
these conditions we may replace the probability by the 
concentration of the reactant at time t, expressed in num- 
ber of moles per m 3 , and describe the mean field chemical 
kinetics of JQ| by reaction-diffusion equations such as, 

d t n A = -k x n A n B (n R + n R *)+k 2 (n R + n R *)+D A V 2 n A 

(2) 



and 



d t n R = kin A n B n R - k 2 n R + D R \7 2 n R 



(3) 



(similar for R*). The rates of reaction, k\ and k 2 have 
units of m 6 s _1 per mole 2 and s _1 respectively. 

These equations can be solved exactly in the 
well-stirred limit (drop spatial derivatives) to obtain 
dn R jdn R * = n R fn R *. This implies that with symmetric 
initial conditions, n R (0) — n R * (0), there can be no chiral 
symmetry breaking, i.e., n R (t) = n R *(t) at all times. A 
fixed point analysis of (J2J and J^J shows that there is a 
line of neutral stability with n A — n R = yfoijki — there 
exist marginal fixed points corresponding to both chiral 
symmetric and asymmetric states. This implies that an 
arbitrary perturbation from the chiral symmetric fixed 
point generated by the slightest stereo-preference, will 
remain as is (i.e., unamplified) and will not give rise to a 
global chiral symmetry breaking @. 

The same analysis may be applied to the reaction- 
diffusion equations if the local concentration 
field is large enough, n A > yfozjkx, there is no sym- 
metry breaking. We next study the effect of includ- 
ing an additive or a multiplicative noise in J5J) and 
Such sources of noise arise in an open chemical system 
through spatiotemporal fluctuations in the number or 



2 



temperature or indeed any other 'fast' variable that cou- 
ples to the chemical rates. The additive (multiplica- 
tive) noise is distributed independently and uniformly 
with zero (nonzero) mean and variance chosen such that 
the concentrations are non-negative at all times. We 
have numerically checked that such an implementation 
of stochasticity does not lead to chiral symmetry break- 
ing starting from symmetric initial conditions. Indeed in 
multiplicative noise models, it is easy to see that the zero 
relative enantiomeric concentration (j) = n R — tir* = is 
an absorbing state 0. 

An important feature of mean-field chemical kinetics 
is that in every interval At, all reactions take place si- 
multaneously. Thus the stochasticity considered so far is 
incapable of generating global chiral symmetry breaking 
from symmetric initial conditions. 

Chemical reactions however have another source of in- 
herent stochasticity, since chemical changes occur during 
molecular encounters which are fundamentally discrete 
stochastic processes; the probability of a reaction occur- 
ing is related to the probability of the reactants encoun- 
tering each other in a time interval At derivable from ki- 
netic theory. This suggests that the appropriate stochas- 
tic description is to treat the reaction-diffusion kinetics 
as a markov process described by a master equation for 
the probability distribution of the number of molecules 
of a given species within a 'local' coarse-grained volume. 
This treatment should reduce to the mean-field model 
when the number of molecules per unit volume is large 
enough to have a well-defined concentration field varying 
smoothly and deterministically with space-time. 

We first construct a local coarse-grained volume by 
noting that at scales smaller than £ ~ y/Dt r , where D is 
the typical diffusion coefficient and t r ~ (kin 2 ) -1 is the 
typical reaction time, the chemical system may be consid- 
ered well-stirred. To incorporate spatial fluctuations over 
scales larger than £, we divide the chemical pool into M 
blocks, each of size £ situated on a cubic (square) lattice. 
The contents in each box is assumed well stirred, however 
each box is connected to its z nearest neighbours allow- 
ing the reactants to diffuse into neighbouring boxes. The 
evolution of the probability distribution of the number of 
molecules of a given species within each box, is specified 
by writing down transition probabilities obeying detailed 
balance, which involve both reaction processes in each 
box and diffusion processes connecting nearest neighbour 
boxes. An efficient stochastic algorithm for solving the 
reaction master equation is the Gillespie algorithm 
which takes as inputs, the initial numbers of each species 
and the chemical rate constants k\ — ci£ 2d and k%. The 
number of molecules of each species is updated whenever 
a reaction (following a collision) occurs; this is a stochas- 
tic event and the probability of its occurence is computed 
from kinetic theory. 

Diffusion can be modeled by transferring a molecule 
of a given species from box j to a neighbouring box j 
with a rate given by A = -D/£ 2 , where D is the diffusion 
coefficient of the given species 0. Thus if N R (t) is the 



instantaneous number of R molecules in block i, then we 
may augment the master equation by a loss term 



P 



Nji (t + At) = N [ R l > (f ) - 1 = AA^ ; (t) At 



r(0/ 



and a gain term 



N R i \t + At) = N%\t) + l =-J2^\t)At. 



We have explicitly checked that this correctly describes 
the diffusive dynamics of R; in the continuum limit this 
recovers the particle current J cx Vn. 

In each box we fix a different random number seed 



(i) 



N 



R* J ' 



The 



and initial data set {N\ 
time evolution of physical observables are computed from 
the stochastic master equation; in all displayed graphs 
we have chosen c\ — 10 _3 s _1 and k% = 10 _1 s _1 and 
have checked that our qualitative results remain unal- 
tered when we vary the rates c\ = 10~ 3 — 1CP 5 and 
&2 = 1CP 1 — 10~ 2 (changing these rates merely changes 
the time and space scales). We have made time t di- 
mensionless by multiplying with k^. Averages denoted 
by (. . .) are over 10 realisations of random number seeds, 
further the temporal behaviour of various observables are 
smoothened by binning over a time interval of t = 5. 

First the diffusionlcss limit, A = : the boxes are un- 
coupled and independent of each other. We follow the 
reaction dynamics in a typical box; the initial data and 
the chemical rates are perfectly chiral symmetric, fur- 
ther recall that the mean field steady state is also chiral 
symmetric. Despite this, in the stochastic evolution the 
populations of R and R* stay roughly the same for a 
while before exhibiting a definite chiral symmetry break- 
ing beyond a time scale r (Fig. 1). Subsequently, the 
population of one of the enantiomers, say R*, goes to 
zero, while the other rapidly takes over the entire pop- 
ulation, (Nr) oc 1 — exp(— (3t), with a rate that is an 
increasing function of ki . The average number of A (and 
B) molecules decrease with time to an asymptote which 
is a weak function of initial conditions. 

Chiral symmetry breaking is best expressed by the 
behaviour of the relative enantiomeric excess $ = 
( Nr+NrI ) w hi° n shows a clear bifurcation before sat- 
urating to ±1 (Fig. 1 (inset)). We compute P(t), the 
distribution of the time scale beyond which chiral sym- 
metry breaking occurs (Fig. 2) : P(t) ~ rexp(— r/r ) 
with a mean (r) that scales linearly with N, for large 
N. Since chiral symmetry is broken in every box, the 
saturation value of the spatially averaged $(£ = oo) — 
A/ -1 J2i ®i(t — oo) = ¥ is zero. 

This dynamical bifurcation may be readily understood 
in a simplified reaction scheme, A ^ R, A ^ R*, 
with the forward autocatalytic reaction rate given by 
kin R /(n A + n R + n R .) and k\nR*/{riA + n R + tir*) re- 
spectively, where it is easy to show the following 10] : 
(i) the mean field equations do not show chiral symmetry 
breaking, (ii) the master equation shows chiral symme- 
try breaking and (iii) the master equation goes over to 
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the mean field equations in the limit of large number of 
reactants. 

Turning on a small A couples the boxes to one another. 
The steady state of this coupled map lattice depends on 
two time scales — the typical chirality breaking time r 
and the diffusion time Td(X). Clearly when Td 3> r, i.e., 
small A, each asymptotic $i is still ±1 with equal proba- 
bility leading to the order parameters, $ = 0, <i> 2 = 1, a 
locally chiral but globally racemic state. 

At intermediate values of A, when r ~ Td, one might 
expect a global chiral symmetry broken state; that this 
is so is seen from Fig. 3 — the magnitude of the or- 
der parameter |<i>| shows a distinct "jump" at A*, which 
gets sharper as M is increased. We have checked that 
this symmetry broken state persists by simply running 
the simulation for longer times. Once the system has 
locked into this symmetry broken state, fluctuations aris- 
ing from diffusive transport can no longer restore the 
symmetry. Note as in the dynamics towards any symme- 
try broken state (e.g., a quench from the paramagnet to 
the Ising ferromagnetic state), an average over different 
random seeds averages over the two degenerate symmetry 
broken configurations. The evolution of the order param- 
eter with time shows a typical nucleation profile (Fig. 4), 
with an initial lag time and a subequent exponential rise 
after a nucleation time r„ which is independent of system 
size (indicating a finite correlation length at the transi- 
tion) and increases with increasing A. The stability of 
this symmetry broken phase can be tested by creating a 
'droplet' of size L containing the opposite phase; we find 
that the mean time for droplet 'evaporation' (r e ) scales 
as L 2 (Fig. 4). 

To understand this novel symmetry breaking transi- 
tion, we coarse-grain our lattice of boxes into cells of size 

the coarse grained <&' = m _1 Y^T=\ ^ where m = £'/£. 
At steady state, takes values 0,±1. In terms of these 
coarse-grained variables, the equations for <&' is given by 

^ = DV' 2 &-V(&)+r,'(x',t) (4) 
at 

where V' = dV/d&' is derivable from a 'fluctuation- 
induced potential' of the form V — r$' 2 — u<f>' 4 + v& 6 , 
with positive coefficients. Discretising the laplacian con- 
verts (0J into a coupled map lattice; if fluctuations induce 
one box to reach a chiral state, say R, then diffusion of R 
from this box to its racemic neighbour will explicitly bias 
the neighbouring box to R. Thus the diffusion D is the 
ordering influence, leading to a discontinuous symmetry 
breaking when large enough. 

At the other extreme of large diffusivity, A — > oo, the 
time scales r 3> Td, i.e., the molecules make rapid ex- 
cursions to other boxes. Thus in the limit of large M, 
the system behaves as a single box with a large number 
N = of molecules. This is the mean-field limit 

where all the reactions happen simultaneously, the fluc- 
tuations leading to chiral symmetry breaking are absent. 
Thus, we expect the steady state to be $ = 0, $ 2 = 0, 



a local (and global) racemic state. We believe that the 
transition from this chiral symmetry broken phase to the 
racemic phase at larger A = A c is sharp and discon- 
tinuous, as can be seen in the simpler reaction scheme, 
A f± R, A ^ R* [T^j. Thus we claim that the steady 
state diagram as a function of A, shows a chiral symmetry 
broken phase at intermediate A and a reentrant racemic 
phase at low A, separated by a sharp discontinuous phase 
boundary. 

The underlying reason for the chiral symmetry break- 
ing in this autocatalytic reaction comes from the stochas- 
ticity in molecular encounters — fluctuations resulting 
in a run of successive i?-reactions, will be followed with 
higher probability by more i?-reactions (autocatalysis). 
Further, because of the autocatalytic nature of the reac- 
tions, once a symmetry broken chiral state has been cho- 
sen, there are no finite-size fluctuations which take the 
system to the other degenerate state. Note that the sym- 
metry broken steady state is not the equilibrium state, 
which can be obtained from the stationary solutions of 
the chemical master equation. The equilibrium distribu- 
tion is clearly chiral symmetric, depends only on the ratio 
of the reaction rates kx/k-z and is independent of the dif- 
fusion coefficients. This fluctuation induced symmetry 
breaking transition occurs in finite size systems and is 
independent of spatial dimensionality. 

To summarize, we have identified a new role for fluctu- 
ations in generating a symmetry broken state in systems 
which are perfectly symmetric in the absence of fluctu- 
ations. We illustrate this phenomenon in the context of 
autocatalytic reaction-diffusion systems : spatiotemporal 
fluctuations drive the chemical system to a chiral sym- 
metry broken steady state. 

We end with some observations — (i) The symmetry 
breaking arising from stochastic switching of the degen- 
erate autocatalytic reactions is a chemical realisation of 
Parrondo's games (ii) Since the diffusion coefficient 
is a function of temperature, our phase diagram allows 
for an interesting interpretation in the context of emer- 
gence of homochirality in the primordial chemical pool. 
At high temperatures, when A is high, the pool is in a 
local racemic state; as the primordial pool cools down, so 
that A < A c , global enantiomeric selection occurs. Once 
selected, homochirality stays on, even when the tempera- 
tures reduce further. Note contrary to popular belief j^|, 
antagonism is not a necessary condition for chiral sym- 
metry breaking, (iii) Is there an analogue of a 'symmetry 
breaking field' by which we may robustly bias the selec- 
tion of a specific enantiomer ? An obvious way is to ex- 
plicitly couple the chemical reactions to an external field 
such as a substrate or circularly polarised radiation which 
would provide a stereospecific bias. Another possibility 
is to drive the chemical soup away from thermodynamic 
equilibrium (e.g., mechanical stirring j0,El), providing 
a chemical ratchet capable of rectifying dynamical fluc- 
tuations to bias 'motion' in chemical space, (iv) The au- 
tocatalytic reactions share features in common with 
the molecular quasi-species model |l4j . namely, replica- 
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FIG. 1: Number of reactants versus time for a given set of initial data : while A and B(o) saturate to a constant population, there 
emerges a definite enantiomeric selection between R(o) and R*(D). Inset shows the time evolution of the relative enantiomeric 
excess {$}, with error bars from different noise realisations. 



tion (R — ► 2R) and mutation (R «-> R*). Fluctuation 
induced selection could be a good candidate for neutral 
mechanisms of selection in the absence of fitness pres- 
sures [la ]. 
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FIG. 2: Scaling of the mean time for chirality discrimination with total number of reactants for different initial data, showing 
roughly (r) oc N. Inset shows the distribution of r has an exponential tail for a fixed initial data, (Na = 100, Nr = 10), 
averaged over 500 realisations to obtain smoothened data. 
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FIG. 3: Saturated value of the magnitude of the spatially averaged enantiomeric excess as a function of A shows a jump at A* 
which gets sharper with increasing M. 
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FIG. 4: Time evolution of (|$|) showing typical nucleation and growth, for various system sizes M and diffusion rate A. Note 
that the nucleation time is independent of M and decreases with increasing A. Inset shows the scaling of the 'evaporation' time 
of a single droplet of size L, (r e ) ~ L z , a best fit gives z m 1.93 ± 0.2. 



